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The scaling behaviour of euclidean quantum gravity at an asymptotically safe critical point is 
studied by means of the exact renormalisation group. Gauge independence is ensured via a specific 
parameterisation of metric fluctuations introduced in a recent paper. Within a non-perturbative 
approximation the beta function for Newton’s constant takes a simple form to all orders in h. A 
UV fixed point is found to exist for d <7 spacetime dimensions at which the critical scaling can be 
assessed. The critical exponent for the Newton’s constant v is found to be regulator independent 
close to two dimensions. Applying Litim’s optimisation criteria we find v rs 1/3 in four spacetime 
dimensions. This value is in agreement with lattice studies supporting the existence of a second 
order phase transition between strongly and weakly coupled phases. 


Introduction .— An important open question in the¬ 
oretical physics is whether a consistent renormalisable 
theory of gravity exists in d = 4 dimensions. Asymp¬ 
totic safety [T. offers one possible realisation of such 
a theory, whereby the ultra-violet (UV) limit is con¬ 
trolled by an interacting fixed point of the renormal¬ 
isation group (RG). The essential idea of asymptotic 
safety is intimately related to critical phenomena in 
systems where a large number of collective degrees of 
freedom play a role. In particular, interacting fixed 
points generally occur where the physical system un¬ 
dergoes a second order phase transition [2]. Thus, as 
in critical phenomena, a fixed point is characterised 
by a set of universal exponents. In quantum gravity 
such an exponent v can be identified with the diver¬ 
gence of an RG invariant correlation length 
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as the (dimensionless) bare Newton’s constant Gb is 
tuned to its critical value G* where A is the UV energy 
scale (with Gb = G(A) = K d ~ 2 G\ where G\ denotes 
the dimensionful bare coupling). Lacking any exper¬ 
imental means to measure such scaling behaviour, it 
is crucial that v can be computed by complementary 
theoretical approaches. 

Early progress was also made by exploiting the 
e = d - 2 expansion whereby at two-loops, and up to 
order e 2 , the value 1/v = e+ |e 2 has been found (3]. 
In four dimensional quantum gravity two main tools 
have been used to study critical behaviour, namely 
lattice models II 0 and the exact renormalisation 
group (ERG) [BUTTl. On the lattice critical scaling 
can be examined via quantities such as the mean cur¬ 


vature 

(R) = (f d^R^/V (2) 

where V is the d-dimensional volume and R is the 
scalar curvature, which is integrated over spacetime, 
with all quantities in units of the lattice spacing 
I = A -1 . One may then use the fixed point scaling 
properties of free energy F{Gb ) = -log Z(Gb)/V to 
infer that JT2] 

F(G b ) ~ C d ~ ]G* - G b \ dv as G b - G*, (3) 

where for gravity Z(G b ) = / T>ge 16 ' ,G b f d 1 \G R+ --- j g 
the partition function. The scaling behaviour of (R) 
can then be determined by differentiating with re¬ 
spect to the bare coupling G b to arrive at (R) = 
l&nGl-^F ~ |G* - Gfcl*'" 1 ~ £ 1/!y " d . By studying 
the scaling behaviour of such quantities, for a lat¬ 
tice model based on Regge’s simplicial formulation of 
gravity the value gua, 


i'Lat » 0.335(4), (4) 

has been found in d = 4 dimensions, where the error is 
due to statistical uncertainties. In addition arguments 
based on the non-local effective field equations m 
imply that only v = 1/3 leads to consistent solutions. 
In this letter we will use the ERG to compute v finding 
a value 


^erg ~1/3. 


( 5 ) 


in d = 4 for suitably optimised regulator schemes. 
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Gauge independence and the one-loop beta 
function. — One issue that must be resolved is the 
possible gauge dependence of z/, which suggests that 
unphysical degrees of freedom contribute to the beta 
functions from which |l]) can be obtained. In a recent 
paper a parameterisation of metric fluctuations was 
put forward for which gauge independence is ensured 
at the one-loop level DSI- The metric 7 ^ is given 
in terms of fluctuating fields around a background 
metric g^, taking the explicit form 

2 

( (j\ d h X ~ O’ 

1 + ”2 I 9^x\ e L / = Qiiv + ft/xi/ + + ••• ? (®) 

where h is a trace-free symmetric matrix with compo¬ 
nents \h\ gl , = and a is a scalar field which param- 
eterises the conformal fluctuations. Gauge indepen¬ 
dence is achieved since the hessian of the bare action, 
expanded around an arbitrary Einstein background, 
does not involve terms proportional to the equations 
of motion. This property follows from © since the 
volume element ^7 oc 1 + ^ is linear in a. One then 
observes the cancelation of gauge fixed fields as well 
as the conformal fluctuations with the measure of the 
one-loop functional integral. This leads to the fol¬ 
lowing beta function for the dimensionless Newton’s 
constant in d > 2 dimensions pb] 

Pg = {d- 2)G - jj(18 ~ N g ) G 2 , Ng =^^-,{ 7) 

universal up to a rescaling of G, where N g denotes the 
polarisations of the graviton. For d = 2 the contribu¬ 
tions of the conformal fluctuations a are absent and 
one finds /3(G) = -|26 G 2 reflecting the 26 dimensions 
of string theory. 

Given the beta function for Newton’s constant the 
value of v can be obtained from /3(G) by differentiat¬ 
ing: 

- d - 2 + quantum corrections ( 8 ) 

G=G* 

where G* is the values of G at a fixed point for which 
/3(G*) = 0. One may then integrate the RG flow in 
the vicinity of G* whereby the correlation length ([I]) 
enters as an integration constant. At one-loop we 
simply reproduce the canonical scaling 1 /zz = d - 2 
which follows from the dimensionality of the (inverse) 
Newton’s constant. 

Flow equation- Here we shall study the flowing ac¬ 
tion Tfc[y>] which depends on the momentum scale k 
and generalises the effective action T to which it is 



equal in the limit k -* 0. An important property of 
Tfc [y?] is that it obeys the exact flow equation [17 [ fl8] 


5 t r fc = 



if } ’ 


( 9 ) 


where t = log(fc/A) denotes the RG time. In the case 
of gravity [5] the fields ip = {h^ v ,a,...} denotes the 
metric fluctuations fields as well ghost fields. The 
right side of flow equation is a super-trace involving 
the hessian of the action I'f and an infra-red (IR) 
regulator IZk- Utilising ([ 9 ]) will allow us to compute 
quantum corrections to the critical exponent ([s]) . 

The details of our calculation proceed along the 
lines of |16l , taking advantage of the parameterisation 
©. but where we now consider the flowing action T k , 
rather than the functional integral from which it can 
be derived [S]. In particular we consider an action T*. 
of the Einstein Hilbert form, 

n = /kVy(A*-^) + ..., do) 


plus gauge fixing, ghost and auxiliary terms from Ja- 
cobians in the functional measure. Here Gk and 
denote the fc-dependent Newton’s coupling and the 
vacuum energy respectively. Exploiting the parame¬ 
terisation ^ we can then find the RHS of ^ upon 
setting all fluctuating Helds to zero. This constitutes 
a background field, or single metric, approximation 
to which we will confine ourselves, noting that meth¬ 
ods that go beyond this [Hfl - E?! could also be applied. 
Compared to the one-loop calculation the background 
field approximation implies an RG improvement lead¬ 
ing to a beta function to all orders in G. 

To perform this RG improvement we assume that 
the hessian of ghosts and auxiliary fields depends 
on Gk such that the unphysical degrees of freedom 
continue to cancel as in the one-loop approximation. 
These cancellation includes the conformal fluctua¬ 
tions and implicitly the non-propagating transverse- 
traceless fluctuations with the Jacobians. The flow 
equation © then takes a form identical to the one- 
loop flow but with the bare Newton’s constant Ga 
replaced by the running coupling Gk 


d t T k = 


1 ®tR-k, 2 lrp r 9tTZk,i 

2 r rg + ft fci2 2 r r g + 


(ii) 


where the first term is a trace over transverse-traceless 
fluctuations hi,, and the second term is a trace over 
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transverse vectors arising from the functional mea¬ 
sure. The hessians appearing in © are given by 

r S = ( 16 ^ G fc)' lA »> 

A 2 <yV, = (-V 2 w -2 R^ a v P<p al3 ) , (12) 

= ( _ V 2 ^ - Rfj, v ) 

with R/tavp and V 2 denoting the Ricci and Rie- 
mann curvatures and the laplacian respectively. The 
regulator 7 Z k takes the form 



rg(A„) + n k , n (A n ) = rg(A„ - p 2 ), (13) 

with Pi E An + k 2 C(AJk 2 ) . 

Here C(z) is a dimensionless regulator function with 
the following properties 

C(l) = 1, C(z -*• oo) -* 0, 0 + G(0) = finite, (14) 


but is otherwise an arbitrary monotonic function. 
The first condition is simply an arbitrary normalisa¬ 
tion condition which removes the degeneracy of C(z ) 
under a rescaling k 2 -*■ ck 2 for which observables are 
invariant. The condition that the regulator vanishes 
for large momentum ensures that it is a Wilsonian 
regulator whereby the high energy modes are unsup¬ 
pressed, while for low momentum the regulator be¬ 
haves like a mass 7Z k ~ k 2 . 

Beta functions .— Using the early time heat kernel 


expansion for the operators (12) one can evaluate the 


traces on the right side of (11). Comparing terms 


up to linear order in the curvature we obtain the beta 
function d t G = (3g = G(d-2+r]G) for the dimensionless 
Newton’s constant G s k d ~ 2 G k with 


V G = 


2 (N g - 18) GI d / 2 -i 


3 (47r) 2 T(| - l) + (N g - 18) GI d / 2 -% 


(15) 


and the beta function 9 ( A = (3\ for the dimensionless 
vacuum energy A = k~ d X kl given by 


Figure 1: The critical exponent v plotted as a function 
of the regulator parameter s which modifies the renormal¬ 
isation scheme. The corresponding inverse propagators 
P 2 (z) are plotted in the inset with values of s ranging 
from s = s m in to s = oo (bottom to top). 


We note that the quantum corrections to /3g continue 
to be proportional to the universal factor N g - 18 
found in the one-loop approximation. 

UV critical scaling .— Solving for /3(G*) = 0 one 
finds a fixed point for positive G* for N g < 18 and 
hence in d < 7 dimensions. At this fixed point the 
critical exponent for Newton’s constant (|8| reads 

l/v=(d-2) + (d-2) 2 ^^. (17) 

^ x d/2-l 


This formula is in fact independent of the universal 
factor | (N g - 18) which appears in /3(G) and arises 
from the heat kernel coefficients of Ai and A 2 . In¬ 
stead © owes its form to the general scaling proper¬ 
ties of the flow equation. These properties are shared 
by a conformally reduced toy model [55] for which the 
critical exponent ( fT7| was also found. We note that 
v is real unlike the critical exponents found for gauge 
dependent RG flows which have been analysed in [ 23 ] ■ 
If we take the limit d -*■ 2 both integrals appearing 
in (17) diverge with as 2/(d-2) and we may therefore 
determine the critical exponent to second order in the 
e = d - 2 expansion 


/3a — ~d\ + 


N g ~ I’ndi) 
(47r)#r(f) 


(16) 


Here the regulator dependent integrals I n and I n take 
the form 


r 

- L "n 



z n ~ 1 C(z) 
+ C(z 



C'(z 
+ C(z) ■ 


1/^ — e + —e 2 + C> (e 3 ) , (18) 

independently of the regulator function G(z). Ex¬ 
trapolating to e = 2 this gives v - 1/4. This result 
can be compared to the 2-loop result [3] which gives 
1/v = e + |e 2 = 4.4 for e = 2. 

Optimised scaling exponent in four dimensions .— 
Away from two dimensions we find that v generally 
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depends on the regulator. In what follows we will 
concentrate on the case d = 4. In this case we have 
computed v for a wide range of regulators G(z), find¬ 
ing values in the range 1/4 < v < 1/2. For example 
the entire range 1/4 < v < 1/2 is found for the one 
parameter family C{z) = a 0(1 - z) (where 0(x) is 
the Heaviside theta function) for which we find 


l/u = 2 + 


2a log (a 1 + l) 


a log (a -1 + 1) + log(a +1) 


(19) 


for 0 < a < oo (although a < 1 is not compatible 
with the normalisation G( 1) = 1). Interestingly this 
corresponds to the entire range of values for which the 
phase transition is second order, as implied by 
with the lower bound agreeing with the e-expansion 


(18). Furthermore we note that the value v = 1/3 is 
attained for a = 1. 

A better determination of u, within the current ap¬ 
proximation, is expected if we apply an optimisation 
criteria to the space of regulators C(z). To this end 
we first note that there must be a non-zero gap in the 
inverse propagator 

H g 2 ap [G] = P\z min ) = G(z min ) + z min * 0, (20) 

where z m - ln is the value of 2 > 0 for which P 2 (z) takes 
its global minimum. Generically we find that it is 
only in the limit where the gap vanishes P 2 ap -*■ 0 
that the value v - 1/4 is attained. To improve our es¬ 
timate of v we will apply Litim’s optimisation criteria 
|25| which is designed to improve the convergence of 
approximate solutions where only a finite number of 
operators (e.g. and R) are retained. Specif¬ 
ically, we demand, in addition to (14), that the gap 


(20) is maximised under the regulator scheme (RS), 
which implies that 


^gap[a p t] = max(P g 2 ap [C]) = 2, 


( 21 ) 


where the specific value max(P 2 ap ) = 2 is subject to 


our normalisation (14). Any regulator G = Copt which 


obeys (21) is said to be optimised. This criteria has 


been applied to the 3 d Ising universality class in [26] 
leading to good convergence of the derivative expan¬ 
sion and a small systematic error in universal quan¬ 
tities. A class of such optimised regulators is given 
by 


Copt (2) = SZ b 


1 


(1 + s) z ° - 1 


where b opt (s) = 


(1 + s ) log(l + s) - s ’ 


( 22 ) 


for values of s (& op t) hr the range s m i n < s < 00 (6 max > 
b 0 pt > 0), with s m ; n = 2.512 (6 max = 1.322), for which 


(21) is satisfied for the global minimum. We note that 
in the limit s -»■ 00 we have G op t(~) = 1 /z. The values 
oil/v obtained for this class of regulators are plotted 
in fig. [I] as a function of s. One observes that 


^opt ~ 1/3 


(23) 


for all permissible values of s with a variation on the 
level of ~ 1%. For s -*■ 00 we obtain v - 1/3 ex¬ 
actly. To check that this is not dependent on this 


class of regulators (22) one may also calculate v for 


Litim’s optimised regulator m which has the form 
(2 - ,s)0(2 - z) for which we again obtain v - 1/3. 
We have also evaluated v with other optimised reg¬ 
ulators and in each case found v a 1/3. We observe 
therefore that, although for generic regulators we find 
values 2 < 1/^ < 4, the application of the optimisation 
criteria (21) predicts values v rj 1/3 which agree quan¬ 


titatively with the lattice result Q. 

Scaling of the free energy .— We now wish to make 
contact with the scaling arguments made in the in¬ 
troduction and in particular recover the scaling law 
for the free energy ([3|. This can be obtained from 
scaling of physical quantities at k = 0 in the limit 
Gb -*■ G* where Gb = h. d ~ 2 G\. Provided Gb < G* 


the beta functions (15) and (16) for G and A define 


a flow which leads to a classical regime at k = 0 for 
which the dimensionful couplings flow to constants 
Gfc^o = Gjv and Afc^o = Ao- These trajectories are 
then well defined for all A > k > 0. An individual tra¬ 
jectory is characterised by the dimensionless number 

d 

To = 8nGf^ 2 Ao obtained in the IR limit. We can then 
obtain the free energy from the effective action T = T 0 


F = r [g^]/V = 


TO 


2An 


d ~ 2 47r(d-2)G^- 2 


(24) 


where F is evaluated on a solution to the equations of 
motion g p „. One may then integrate /3(G) between 
G = Gb = and G = 0 to find the following scaling 
relation 


Gjv = A 2 d e exp 


r Gb 1 
(d- 2) J dx 


P(x). 


£->0,G b ->-G* 


5 

= A 2 - d Gj G *~ Gb ) 

b = b 0 pt 

\ G* ) 


-v(d- 2 ) 


+ ... ~ 


d -2 


(25) 


Inserting this into (24) we recover the scaling law 
([3]) . In addition this scaling indicates that the Planck 
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length Ip] e G ^ -2 oc £ emerges due to the existence of 
the fundamental length scale £. 

Discussion .—In the lattice theory [3] for which 0 
was obtained the weak coupling phase Gb < G* gives 
rise to breached polymer phase which is usually un¬ 
derstood as an effect of the conformal instability of 
the Euclidean action. Here we instead find a well be¬ 
haved weak coupling phase with classical scaling. In¬ 
deed, we have implicitly Wick rotated the conformal 
factor EE] resulting in the cancelation of conformal 
fluctuations with the functional measure. In turn on 
the lattice there is a strong coupling phase G& > G* 
with spacetimes of negative curvature [3 [2?5] while for 
the beta function found here we run into a pole of (15) 
for Gb > G*. We observe therefore that currently the 
euclidean lattice theory only has access to a physical 
strong coupling phase whereas the ERG has access 
only to a weak phase. Since the conformal instability 
is due to the Euclidean nature of the path integral one 
may hope to get access to both phases via a suitable 
Wick rotation of the lattice theory. Such a rotation is 


offered by the causal dynamical triangulation (CDT) 
lattice theory [5], where a phase with classical scal¬ 
ing is observed [31]]. In turn one might hope that a 
strong coupling regime is accessible to the ERG once 
we go beyond the approximation used here e.g. by 
including higher order curvature invariants [91411] or 
utilising vertex expansions USM]- 

We conclude that the quantitative agreement be¬ 
tween lattice and the continuum studies pressed here 
suggests that a second order phase transition exists 
between strongly and weakly coupled quantum grav¬ 
ity. This critical point could then provide a contin¬ 
uum limit for quantum gravity. 


Acknowledgements 

The author would like to thank Jan Pawlowski and 
Daniel Litim for useful comments and Herbert Ham- 
ber for correspondence. 


[1] S. Weinberg (1979), in Hawking, S.W., Israel, W.: 
General Relativity. 

[2] K. G. Wilson and J. B. Kogut, Phys. Rept. 12, 75 
(1974). 

[3] T. Aida and Y. Kitazawa, Nucl.Phys. B491, 427 
(1997), hep-th/9609077. 

[4] H. W. Hamber, Phys. Rev. D61, 124008 (2000), hep- 
th/9912246. 

[5] J. Ambjorn, A. Goerlich, J. Jurkiewicz, and R. Loll 
(2013), 1302.2173. 

[6] M. Reuter, Phys. Rev. D57, 971 (1998), hep- 
th/9605030. 

[7] W. Souma, Prog. Theor. Phys. 102, 181 (1999), hep- 
th/9907027. 

[8] D. F. Litim, Phys.Rev.Lett. 92, 201301 (2004), hep- 
th/0312114. 

[9] A. Codello, R. Percacci, and C. Rahmede, Annals 
Phys. 324, 414 (2009), 0805.2909. 

[10] D. Benedetti, P. F. Machado, and F. Saueressig, Nucl. 
Phys. B824, 168 (2010), 0902.4630. 

[11] K. Falls, D. Litim, K. Nikolakopoulos, and 
C. Rahmede (2013), hep-th/1301.4191. 

[12] F. J. Wegner, Phys.Rev. B5, 4529 (1972). 

[13] T. Regge, Nuovo Cim. 19, 558 (1961). 

[14] H. Hamber, to be published (2015). 

[15] H. W. Hamber and R. M. Williams, Phys.Lett. B643, 
228 (2006), gr-qc/0607131. 

[16] K. Falls (2015), 1501.05331. 


[17] C. Wetterich, Phys. Lett. B301, 90 (1993). 

[18] T. R. Morris, Int.J.Mod.Phys. A9, 2411 (1994), hep- 
ph/9308265. 

[19] D. Becker and M. Reuter, Annals Phys. 350, 225 
(2014), 1404.4537. 

[20] A. Codello, G. D’Odorico, and C. Pagani, Phys.Rev. 
D89, 081701 (2014), 1304.4777. 

[21] N. Christiansen, B. Knorr, J. M. Pawlowski, and 
A. Rodigast (2014), 1403.1232. 

[22] I. Donkin and J. M. Pawlowski (2012), 1203.4207. 

[23] K. Falls (2014), 1408.0276. 

[24] S. Nagy, B. Fazekas, L. Juhasz, and K. Sailer, 
Phys.Rev. D88, 116010 (2013), 1307.0765. 

[25] D. F. Litim, Phys.Lett. B486, 92 (2000), hep- 
th/0005245. 

[26] D. F. Litim and D. Zappala, Phys.Rev. D83, 085009 
(2011), 1009.1948. 

[27] D. F. Litim, Phys.Rev. D64, 105007 (2001), hep- 
th/0103195. 

[28] P. O. Mazur and E. Mottola, Nucl.Phys. B341, 187 
(1990). 

[29] W. Beirl, E. Gerstenmayer, H. Markum, and 
J. Riedler, Phys.Rev. D49, 5231 (1994), hep- 
lat/9402002. 

[30] J. Ambjorn, J. Jurkiewicz, and R. Loll, 
Phys.Rev.Lett. 93, 131301 (2004), hep-th/0404156. 




